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Abstract 

We have recently shown that the large direct CP violation observed in the decay 5' 







^^ \ TT+TT", and the enhanced branching ratio for B^ -^ tt'^tt^, can be understood by invoking 

a small mixing of the tttt system with the dominant DD channel. We examine here the 

analogous effect of DD ^^ KK mixing on the rare decays i?° — > K'^K^, B~ — > K^K^ and 

f^ ' i3° — > K'^K^ . We find (a) significant values for the asymmetry parameters C and S in 

B^ / B^ -^ K^K'^, (b) a possible enhancement of the suppressed mode B^ -^ K^K^, (c) 

f^ ■ a correlation between the three decay channels following from a triangle relation between 

I I amplitudes Aj^f,-^ — A^+x- = A^^-xo- The pattern of asymmetries and branching ratios 

jr~^' is compared with that derived from the short-distance QCD penguin interaction. 

(— I \ 

The rare decays B^ — > K^K^, B~ — > K^K^ and B^ — > K^K^ have a rather simple descrip- 
tion in the conventional approach, based on the short-distance Hamiltonian for b — > dss and the 
^ ' assumption of factorization of hadronic matrix elements ^ [21 ^ . The dominant interaction 

is the QCD penguin transition h — > dg* -^ dss, which leads to the following amplitudes ^ (the 
absolute square of an amplitude gives the branching ratio): 

A{W^K^K^) = iVo[Ai + APAJ 

A{B- -^ K-K^) = A{W -^ K^K^) (1) 

A(B^ ^K+R-) = 



Here At = VtbV*^, Ac = VcbV*^ are CKM factors, and AP is a term of order (tti^ — m'^)/m1 
that we discuss below. The normalization Nn contains the penguin factor ^^In— 5-, as well as 
decay constants and wave-function overlaps. A typical estimate is A^o ~ 0.12, which implies a 
branching ratio By{B^ /W -^ K^K^) « 1.0 x 10"^. The notable feature of Eqs. ((U is the equality 
of the decay rates of B*^ — s- K^K^ and B~ -^ K~K^, as well as equal CP-violating asymmetries 
between W -^ K^W and 5° -^ A'°!^,_and between B' -^ K~ K^ and B+ -^ K+K^. Also 
notable is the suppression of the decay B^ — >■ K^K^ (which persists even after inclusion of 
annihilation contributions, the branching ratio being ~ 10^^ ^). 

* E-mail address: sehgal@physik.rwth-aachen.de (L. M. Sehgal) 



Of particular interest to us is the term AP in Eq. (^, which represents the residual contri- 
bution of cc and uu intermediate states in the QCD penguin diagrams |51ini. This complex con- 
tribution depends on the effective (time-like) q'^ of the virtual gluon that mediates the B -^ KK 
transition. An explicit expression for AP(gr^) is given in Ref. ^, and suggests that |AP| can 
be as large as 0.3 for {q'^)es ~ g^^b- ^^ the limit in which AP is neglected, the asymmetry 
parameters C and S for the decay B^/B^ -^ K^K^, defined as 

^ - 1,|\ _|2 ' -^ - 1,|\ _|2' i^-' 

^°-^° PA{B0^K0K<^) ' p 

vanish identically: C = S = 0. Non-zero values of C and S are thus associated with effects 
of cc intermediate states, embedded in the function AP(g^). The behaviour of C and S as 
a function of q'^ is shown in Fig. where we have used A^ = |Au|e~*''' = 3.6 x 10~^e~'^'' , 
Ac = -|Ac| = -8.8 X 10-3, At = -(A„ + AJ and 2/3 ^ 45°. 

In the present note, we re-examine the amplitudes for the B -^ KK channels, from the 
standpoint of inelastic final-state interactions coupling the KK and DD systems, and suggest 
an alternative to the term APAc in Eq. (^ . The motivation comes from a recent analysis of the 
decay channels B -^ tttt [Jj , in which we showed that a mixing of the strong-interaction channels 
vrvr and DD in the isospin / = state helps to resolve the puzzling observations of large direct CP 
violation in B^/B^ -^ ir^n" |Hj and an enhanced branching ratio for B^ /B^ —^ tt^-k^ [2]. There 
is a natural explanation for these correlated effects in terms of final-state interactions in systems 
of physical hadrons, which respect the strong-interaction conservation laws. Experiment |10j 
shows that DD systems are produced with large branching ratio (i.e. large decay amplitude, 
proportional to Ac). Production of such systems in a real intermediate state, followed by a 
small mixing into the final state, results in an important imaginary contribution to the decay 
amplitude. This gives rise to asymmetries, and can enhance an (initially suppressed) branching 
ratio. The same physical idea will be implemented here, via the mixing of DD and KK states 
with the same strong-interaction quantum numbers. Introduction of the DD as an intermediate 
state in the decay B -^ KK acts as a long-distance substitute for the cc contribution in the 
QCD penguin loop parametrized by AP{q'^). Since the discussion is on the level of physical 
hadronic channels, no reference is made to the uncertain variable q'^, which appears in the short- 
distance contribution. Here we deal with the contributions to individual exclusive decays, of 
empirically-established ^Oj) real intermediate states of physical hadrons. We obtain different, 
correlated results for asymmetries and branching ratios. 

Mixing occurs between the following pairs of states, possessing definite strong-interaction 
quantum numbers. 

\kk)+ = ^(k'^'k^ + k+k-)] 

\DD)+^^{D+D-+D^W) j ^-U, C.-+, J -0 (3j 

\kk)^ = ^(k^'k^-k+k-)] 

\DD)^ = ^{D+D- -D^D^) ] ^-^>t.- ,J -U 



As a consequence, the physical decay amphtudes of B^ to these states are |7j 
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SI 
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(4) 



Here At^—, At,—, A^^—, A^— denote the bare decay amphtudes, m the absence of final-state 

KK^ DD^ KK^ DD -^ ^ ' 

1 1 

interactions. The mixing effects are contained in the rescattering matrices S^ and St-, where 
5+ and S- are the strong-interaction S matrices between the pairs of states given in Eqs. (jSJ- 
In the limit in which the elastic phase shifts for {KK)± -^ {KK)± and {DD)± — > {DD)± are 

neglected, the rescattering matrices S^ reduce to the simple form ^J 



si 



cosb 
isin( 



zsm( 
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(5) 



where 0+ and 9- are two phenomenological mixing angles. The physical decay amplitudes of 
B^ to the states K^K^ and K^K^ can then be written as 
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cos 9+{A^o-]^ + Ak+K' ) + cos 6*- (i^o;^ - ^r+r- ) 
sin9+{Aj:,+D- + A^^q-^) + s\d.9^{Aci+D' - ^d^W 

cos 9+{A^o-j^ + Ak+k- ) - cos 6'_ (i^o;^ - ^r+r- ) 
sm 9 +{Aj:,+j:,- + i^o;^?)) - sin6i_(i£,+£,- - i^o;DO 



(6) 



For the charged decay mode B -^ K K^ , involving the 1 = 1 state K K^ , mixing with the 
state D^D^, with mixing angle 0_, implies 



^R-RO = cos9^Aj^-j^o +ism9-Aj^~j^o 



(7) 



For the bare amplitudes, we use the leading term of the QCD penguin model expressed by Eqs. 



A 

^KORO 



^R-RO = NoXt, Aj^+j^- = 0. 



(8) 



The bare amplitudes for B^ -^ D^D , B^ — > D^D^ and B ^ D D^ are taken to be those in 
the Bauer-Stech-Wirbel factorization model [21 El 



Ad+ 



D- 



Ad-do = N'Xc, A^ojTo = 



(9) 



where A^' = 1.79 is determined jjj from the measured branching ratio Bv{B^ / B^ —f D^D^) 
2.5 X 10~^ ^ni- With this input, the physical decay amplitudes, for small mixing angles (cos 9± 
1, sin 0-1- ~ 9±) take the form 



A ~ 

^K'iRO - 



A 



R+R- 



NoXt + iN'X,-{9+ + 9^) 
iN'X,l{e+-9^) 



(10) 
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The amplitudes satisfy the triangle relation (for all 6±) 

'^ko'ko ~ ^K+K- = ^K-RO (11) 



In what follows, we discuss the observable consequences of the amplitudes in Eqs. (|10|) and 
compare them with those obtained from the short-distance penguin model expressed by Eqs. 

1 Branching Ratios (averaged over B and B) 

From Eqs. (|1U|) . we obtain the average branching ratios 

Bt{b^/W^k-^k~) = N'^\x,Ae+-e^f (12) 

Br{B-/B+ -^ K~K^/K+W) = |iVoAtP + N''^\\c\^e'i 

Notice that for different values of ^+ and 0_, the branching ratios of B^ /B^ —> K^K^ and 
B~ /B^ —>■ K~ K^/K~^K^ are no longer equal. Futhermore, for 6^ ^ 6-, the decay rate of 
B^/B^ — > K~^K~ does not vanish, in contrast to the expectations from Eq. (^. The correlation 
of the three branching ratios is shown in Fig. 121 which gives the relation between the K^K^ 
and K^ K^ / K~^ K^ rates for given values of the K~^K~ branching ratio. The empirical limits on 
these decays are Bi{B^ /W -^ K^W) < 1.6 x 10"^ [l2], Bi:{B°/W -^ K+K') < 0.6 x 10"^ ^^ 
and Br{B-/B+ -^ K-R^/K+K'^) < 2.2 x 10"^ IT31. These data constrain the mixing angles 
d± to the intervals 16*+ + 6l_| < 0.15, \9+ - 9-1 < 0.10 and |6'_| < 0.07. 



2 C and S Parameters for B^/B^ -^ K^K^ 

The asymmetry parameters C and S for direct and indirect CP violation respectively, calculated 
from the amplitudes in Eqs. (|1UI) . are shown in Fig. 13 For values of |0+ + ^-| of order 0.1, 
the predicted values are \C\ ~ 0.4, S ~ —0.3. To facilitate comparison with the QCD penguin 
model, we show in Fig. ^the trajectory of allowed values in the C — S plane. The model with 
DD <-^ KK mixing produces the trajectory drawn as a solid line, with points labelled by the 
value oi9 = 6j^ + 0^. The penguin model predicts the trajectory drawn as a dashed line, labelled 
by values of g^, the virtual gluon mass. For the theoretically favoured region ^ < :^ < 2 0) the 

b 

short-distance model predicts C ~ 0.2, S ~ —0.2. While the two trajectories intersect, there is a 
large domain in which the predictions of the two models are quite different. Since Aj^+j^- arises 
from the final-state interaction for 9-^- 7^ ^-, we have Cj^+x- = and S^+k- — ~ sin 2/3 = —0.7. 

3 Effects of a Small Elastic Phase in the Rescattering Matrix 

1 

The matrices S^ that connect the bare and the physical amplitudes in Eq. @ can have a more 
general form than that given in Eq. , if the S matrix is allowed to have non-zero elastic phases 



in the KK and DD channels 7^. For small phases, the net effect is to introduce a further phase 
factor e on the second term in the B'^ -^ K^K^ amplitude in Eq. (|10() . The effect of such a 
modification on the trajectories in the C — S plane is shown in Fig. [51 for two choices 5 = ±20°. 
The effect on the average branching ratios is exhibited in Fig. |H1 

4 Comments on other Models for B -^ KK 

Predictions for B -^ KK decays appear in various analyses dealing with B — > PP transitions, 
where P denotes a pseudoscalar meson. In Ref. ^S], an SU{2>) symmetry is assumed, and decay 
amplitudes are written as linear combinations of phenomenological parameters associated with 
different quark topologies. These parameters are estimated from a fit to all B -^ PP channels. 
The model gives equal rates for K^K^ and K~ K^, and the K^K^ rate is strongly suppressed, 
exactly as in the QCD penguin amplitudes given in Eqs. (^. A similar result occurs in the 
QCD factorization approach of Ref. 0|, which estimates an asymmetry parameter for K^K^ 
of C ~ 0.16. In Ref. J16j . the C and S parameters for B^/B^ -^ K^K^ are estimated to be 
C ~ 0.12, S ~ —0.16, similar to Ref. P, using a model like that in Eqs. (^. In Ref. |17j . 
data on i? — > PP decays are fitted using elastic scattering effects. The fits have very large 
elastic scattering phases. Then, they give a branching ratio for the K^K~ channel of up to 
0.6 X 10~^, a large asymmetry parameter C for the K^K^ channel, and also a Cj^+k- of the same 
sign. Smith ^H] has carried out an analysis of i? — > PP decays which discusses SU{3) elastic 
scattering, and separately, effects of PP <-> DD mixing using a parametrization different from 
ours. In the limit of retaining only the inelastic mixing effect, his results for B^ / B^ —^ K^K^ 
resemble ours, but the branching ratio B^/B^ -^ K'^K~ vanishes. In our approach, an enhanced 
rate for K^K~ occurs as long as the mixing angles 6+ and 9- for states of different isopin, G- 
parity, are different. Finally, there are papers discussing "charming penguins", which attempt to 
parameterize charm-related, long-distance aspects of the penguin diagram, but make no reference 
to strong-interaction eigenstates of physical hadrons in real intermediate states. Results may be 
found in Ref. ^^1 and the literature cited therein. 

In conclusion, we note that current experiments are steadily approaching the level of sensi- 
tivity at which the decays B -^ KK should soon be measurable. One awaits with interest the 
pattern of branching ratios and asymmetries that these measurements will reveal. Our results 
here show that this pattern can elucidate further important aspects of the underlying dynamics. 
The unexpected, recent oberservations of enhanced rates, for vr'^vr'^ 9 , ir^p^ |20|l21j. and rjco [22j 
are striking consequences of this dynamics jTl |2S1 • 
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Figure 1: C and S parameters depending on the virtual gluon mass q^ in the penguin model, for 
the decay B^/W -^ K^K^ [T]. 
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Figure 2: Br (5-/5+ -^ K-K^/K+W)_as a function of Br(5°/50 ^ K^W). The parameters 
O.x describe the branching ratio of B^ /B^ — s- K^K~ in units of 10^^. 
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Figure 3: C and S parameters as functions of the mixing angle 9 
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Figure 4: C as a function of S in the decay B^ / B^ -^ K^K^ . Sohd hne: mixing model for 
parameters —0.15 < 9 = 9+ + 9- < 0.15, dashed line: penguin model for 0.1 m? < q^ < 0.9 m?. 
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Figure 5: C as a function of S in the decay B^/B^ -^ K^K^. Solid line: mixing model for 
parameters —0.15 < 9 = 9^ + 9- < 0.15 and two different 6, dashed line: penguin model for 
0.1 mg <g2 <0.9mg. 
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Figure 6: Br{B-/B+ -^ K-R^/K+W) as a function of Bt{B^ /'W -^ K^) for a dynamical 
phase 6 = ±20°. The parameters 0.x describe the branching ratio of B^/B^ — > K~^K~ in units 
of 10^6^ 



